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Abstract 

The Sand Pile Model (SPM) and its generalization, the Ice Pile Model 
(IPM), originate from physics and have various applications in the descrip¬ 
tion of the evolution of granular systems. In this article, we deal with the 
enumeration and the exhaustive generation of the accessible configuration 
of the system. Our work is based on a new recursive decomposition the¬ 
orem for SPM configurations using the notion of staircase bases. Based 
on this theorem, we provide a recursive formula for the enumeration of 
SPM(n) and a constant amortized time (CAT) algorithm for the gener¬ 
ation of all SPM(n) configurations. The extension of the same approach 
to the Ice Pile Model is also discussed. 


1 Introduction 

The Sand Pile Model (SPM) is a discrete dynamic model inspired by real world 
physics problems, namely the dynamic of piles of granular materials such as sand 
or cereals in silos. A first discrete dynamic system formulation of sand piles was 
given by statistical physicists. The Sand Pile Model we consider here is a one¬ 
dimensional simplified model, where sand grains are stacked on a number of 
adjacent columns. 

We now give a precise definition of SPM. A partition is an infinite non¬ 
increasing sequence (sj)i>o of natural numbers with finite support. We denote 
by Part(n) the set of partitions of n, i.e. partitions (si)*>o with ]T\> 0 Si = n. 
We notice that we index the components of our sequences starting with zero. 
A configuration is simply a partition. We define the set of SPM configurations 
with n grains, denoted by SPM(n), to be the set of configurations reachable 
from the initial configuration (n,0,...) with the following evolution rule called 
the FALL rule: 

s = (so, • ■ ■, Si, si +1 , ...)—>• s' = (soj .. •, S; — 1, Sj+i + 1,...) 

whenever s; > sj+i + 2. Clearly, SPM(n) is a subset of Part(n). 

The Ice Pile Model (IPM) is an extension of SPM with the following addi¬ 
tional rule: 

(s 0 , ■ ■ ■ , Si, ■ ■ • , SZ+fc', ■ ■ ■) -A (So, ■ ■ ■ , Si - l,Sz+l, . . • , Si + k',Si + k'+l + !,•■■) 
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whenever si — 1 = sj+i = ... = si+k 1 = Si+fc'+i + 1 for k' < k. This rule, 
parametrized by the integer k. is called the SLIDEfc rule. We define the set 
of IPMfe configurations with n grains, denoted by IPMfc(n), to be the set of 
configurations reachable from the initial configuration (n, 0,...) by applications 
of the rules FALL and SLIDEfc. We also have that IPMfc(n) is a subset of 
Part(n). 

Some results in counting configurations in SPM(n) and IPMfc(n) are already 
known. For instance, in 15], recursive formulae for |SPM(n)| are given based 
on an inductive lattice structure of SPM(n), and in [Ti the generating function 
of |IPMfc(n)| is studied and its asymptotic behavior is given. On the front of 
exhaustive generation, recent results (e.g. [61 E]) provide efficient exhaustive 
generation algorithm, more precisely in constant amortized time (CAT), using 
the dynamics in the evolution of sandpile configurations. However, these two 
lines of research rely on different aspects of SPM, unlike other combinatorial 
structures, whose counting and exhaustive generation are often the two sides of 
the same coin. In this article, we would like to provide a general framework in a 
combinatorial perspective for both counting, and efficient exhaustive generation 
of SPM/IPM configurations. 

In this article, we study a recursive structure of SPM configurations that 
determines a recursive decomposition for them. This recursive decomposition 
is the key to our results on a new recursive formula for |SPM(n)| and a new 
CAT algorithm to enumerate all configurations in SPM(n). In Section 2, the 
basic notion of “staircase basis” for SPM is introduced, and with this notion 
we characterize the aforementioned combinatorial recursive structure of acces¬ 
sible configurations in SPM(n), and we obtain in particular a recursive for¬ 
mula to determine their number. We also present a natural algorithm for the 
exhaustive enumeration of SPM(n) using the recursive structure of accessible 
configurations, and we prove it to be CAT. Our algorithm differs from the 
Massazza-Radicioni algorithm because we only use combinatorial properties of 
the accessible configurations and not properties related to the dynamic of the 
system. In Section 3, we give an intuitive presentation of how our idea can be 
generalized to IPM. We conclude with some discussions of possible directions 
of future work. 


2 Staircase bases and recursive structure 

2.1 Staircase bases 

Let Lb{tl) be the lattice obtained by equipping Part(n), the set of partitions, 
with the “dominance” order ^ defined as follows: for s,i€ Part(n), s -< t <=>- 
Vj, £i=o s* > ELo U- From EJ we know that SPM(n) is a sublattice of Ls{n). 
We can also define the following partial order called sequence order or covering 
order on Part(n): for s,t £ Part(n ), we write s < t if and only if for all i £ N, 
we have Si < L- This order can be readily generalized to the set of arbitrary 
sequences of integers. 

Now we will introduce our notion of staircase basis for SPM. 

Definition 2.1. Staircase bases for SPM For k £ N, we define the staircase 
of order k by s(k) = (k, k — 1,..., 2,1,0,...) £ Lb- More precisely, Vi £ 
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N, s(k)i = max(0, k — i). We define B = {s(k) \ k £ N} to be the set of staircase 
bases for SPM. 

We have the following property that follows immediately from the definition 
of FALL. 

Proposition 2.1. For s £ B, s is a fixed point for SPM, that is, we cannot 
apply FALL on any column of s. 

We now define a parameter of SPM configurations, called staircase width , 
which will be crucial in the following. 

Definition 2.2. For some t £ SPM(n), we define its staircase width as the 
integer sw(t) = ma x s /f.)<t k. Furthermore, if sw(t) = w, we call s(w) the 
staircase socle (or simply socle) of t. 



Figure 1: t = (6, 6, 3, 3,1,1, 0) with sw(t) = 5 and with its socle in white 

The staircase width of an SPM configuration is “monotone” with respect 
to the evolution rule FALL, as showed by the following theorem which relies 
esssentially on Proposition 12.II 

Proposition 2.2. For a,b £ SPM(n ) such that a b, we have sw(a) < sw(b). 
More generally, if a A b, then sw(a ) < sw(b). 

Proof. We deal with the case a —> b, the general case follows as a consequence. 
Suppose we had sw(a) > sw(b). In this case, there exists some index i such that 
> s(sw(a))i but bi < s(sw(a))i. This is possible only when we apply FALL on 
a at index i, thus at = bi + 1 and a* = s(sw(a))i. However, cq+i > s(sw(a))i+ 1 . 
By Proposition ^.II FALL cannot be applied at index i, which is a contradiction. 
Therefore sw(a ) < sw(b). □ 

In 2]> th e following characterization of the (unique) fixed point in SPM(?r) 
is given. 

Proposition 2.3. In SPM(n), the unique fixed point with respect to rule FALL 
is 

<f(n) = (k, k — 1 , ..., l + 1 , l, l, l — 1 ,..., 2 , 1 , 0, .. .) 

where (k, l ) is the unique pair such that 0 < l < k and n = ^k(k + 1) + 1. 

As a corollary of Propositions above, we obtain an upper bound for the 
staircase width of elements in SPM(n). 

Proposition 2.4. For all t £ SPM(n), we have sw(t) < \[2n. 

Proof. By Proposition 12.21 and 12.31 sw(t) < 8w{<j>(n)) = k < v / 2 n. □ 
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We can partition the set SPM(n) according to the staircase width of con¬ 
figurations. We define SPM(n,w) = {s £ SPM(n)|su;(s) = w} be the subset of 
SPM(n) of all elements with staircase width w. For w running from 1 to > 

all SPM(n,ui) partition SPM(n). From now on, we concentrate on SPM(n,u>) 
instead of SPM(n) as a whole. To generate SPM(?r) exhaustively, it suffices to 
provide a CAT algorithm to generate elements of SPM(n,u;), with the param¬ 
eter w varying from minimal value 1 to maximal value less than \/2 n. 

2.2 Recursive structure 

In [3j, the following characterization of elements in SPM(n) is given. 

Theorem 2.1. A partition s is in SPM(n) if and only if none of the following 
patterns (also called forbidden patterns) occur: 

• p,p,p for p > 0 (that is, three columns containing the same number of 
grains) 

• p,p,p — l,p — 2,..., q + 2, q + 1, q, q for p > q > 0 (that is, two plateaux, 
one of height p and one of height q, separated by a perfect staircase) 

Using this theorem, we can bound the number of non-zero components of 
an element in SPM(n, u>). Given two finite sequences a,b , we denote their 
concatenation by a ■ b. 

Lemma 2.1. The largest (in the sequence order <) SPM configuration s with 
sq < w is w ■ s(w ). 

Proof. Let s be an SPM configuration with so < w. Theorem 12.11 implies that, 
for any integers z > 0 and k > 0, Si+k < Si — k + 1, thus Sk < w — k + 1 
and s < w ■ s(w). We conclude by noticing that w • s(w) is also an SPM 
configuration. □ 

Proposition 2.5. For a £ SPM{n,w), we have a w+ 1 = 0. 

Proof. For a £ SPM(n, w), we have sw(a ) = w, thus there exists an index i < w 
such that dj = w — i. The suffix (a*, a,+i,...) is also an SPM configuration. By 
Lemma 12.11 we have a^+i <0. □ 

As a consequence, we can express an element a £ SPM(n, w) as a (w + 
l)-tuple a = (ao, • • ■ ,a w ). We now introduce a representation of elements in 
SPM(n,w) obtained by “removing” the socle of each element. 

Definition 2.3. Reduced form for SPM(n, w) For s £ SPM(n, w), we define 
red w (s) as follows: 


Vz £ {0,..., w}, ( red w {s))i = s t - w + i. 


Thus red w (s) is obtained by simply substracting (w, w— 1,..., 0) pointwise from 
s. We call red w (s) the reduced form of s. 

Note that, in the above definition, the subscript w is not necessary. It has 
been added for emphasizing the staircase width of the original configuration, as 
well as the number of components of the corresponding reduced form. Compo¬ 
nents in a reduced form may be zero. We notice that the leftmost component 
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that is equal to zero may be the component in position w. Think for instance 
of configuration (4,3) whose reduced form is (2,2,0). For our proposes, how¬ 
ever, it is important to highlight this apparently superflous component, as the 
leftmost zero component of a reduced form is pivotal to define our recursive 
decomposition. 

Let R(n,w ) be the set of reduced forms of elements in SPM(n, w). The 
map red w is clearly a bijection between R(n,w) and SPM(n,w). We should 
also notice that a reduced form is not necessarily an SPM configuration. An 
example of an SPM configuration and of its reduced form is given in Figure [2] 
where the white part is the socle and the gray part corresponds to the reduced 
form. 


t = (6,6, 3,3,1,1, 0) with 
r = red 5 (t) = (1, 2,0,1,0,1) 

Figure 2: An example of SPM configuration and of its reduced form 

The characterization of elements in SPM(n, w ) translates into a characteri¬ 
zation of their reduced forms. 

Proposition 2.6. A (w + l)-tuple r = (ro, r ±,..., r w ) of natural numbers is in 
R(n , w ) if and only if the following conditions are satisfied. 

(i) £r=o n = n- \w{w + 1). 

(ii) There exists at least an index ig with 0 < io < w such that ri 0 =0. 

(Hi) For all i,j such that 0 < i < j < w, we have r* > Tj — 1. 

Proof. Let t £ SPM(n,w) and r = red w (t). Conditions (i) and (ii) follow 
directly from the definition of SPM(n,u>) and of red w . For condition (iii), we 
notice that Theorem 12.II implies that ti + k < ti — k + 1 for any k > 0 and * > 0, 
thus for any i , j with 0 < i < j < w, rj+w—j = tj < ti — (j—i)+ 1 = ri+w—j+1, 
and we have rj < r» + 1. 

For the other direction, let r be a [w + l)-tuple that satisfies (i), (ii) and 

(iii), and let t be the (w + l)-tuple obtained by adding r and s(w) component¬ 
wise. To prove that r £ R(n,w), it suffices to prove that t £ SPM(n,w). From 
condition (iii), it follows that ti = ri + w — i > ri+i — 1 + w — i = ij+i, thus t is 
a partition. By Theorem 12.II we only need to prove that no forbidden pattern 
exists in t. We suppose that such a pattern exists between column i and j, we 
then have U — (j — i) + 1 < tj, thus ri + w — j + 1 < rj + w — j, which implies 
Ti < rj — 1, contradicting condition (iii). Hence t £ SPM(n,w). □ 

We notice that condition (i) only ensures the weight n to be correct. There¬ 
fore, any ( w + l)-tuple r verifying (ii) and (iii) must be in some R(n,w) for an 
appropriate n. We also notice that condition (iii) is preserved by taking pre¬ 
fixes. That is to say, if a tuple r verifies (iii), then all of its prefixes also verify 
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(iii). With these remarks, we provide the following decomposition theorem for 
reduced forms, which is the main result of this article. Both our enumeration 
formula and our exhaustive generation algorithm rely on this theorem. 

Theorem 2.2. Decomposition of reduced forms 

A reduced form r £ R(n , w) can be uniquely decomposed into the following 
form: 

(r 0 , ■ ■ ■, ruf) — (to, • - *; ti— i, 0, uq, ■ • ■, u w —i —\), 

such that t = (to, • ■ ■, tz-i) and u = ( uq ,..., u w -i- 1 ) verify the following condi¬ 
tions. 

• If u is not empty, for all i £ {0,... ,w — l — 1}, we have Ui £ {0,1}. 

• If t is not empty, for all i £ {0,...,/ — 1}, we have ti > 0. 

• In the case that t is not empty, let m be the minimum ofti for 0 < i < l— 1 

and r' = (to — m, t\ — m ,... ,t;_i — to). Then the tuple r' is in some 

R(n',1 — 1), where n! can be easily calculated. 

We refer to this decomposition by writing r = (l,u,m),r'. We represent an 
empty tuple by a pair of parentheses (). When t is empty, we take to = 0. 

Proof. We start by the validity of our decomposition. As r £ R(n,w ), by 
condition (ii) in Proposition 12.61 there exists some index i such that r* = 0. 
To ensure all parts of f to be strictly positive, we take l to be the minimum to 
such that ri 0 =0. By condition (iii) in Proposition 12.61 we deduce that u is a 
sequence of 0’s and l’s. The tuple r' clearly has only positive parts and has at 
least one 0 component, thus it verifies (ii). The fact that r' satisfies condition 
(iii) is provided by the fact that property (iii) is invariant not only by taking 
prefixes (so t satisfies it), but also when a constant is substracted from every 
part of a partition (and r' is obtained from t by subtracting the integer to from 
all its parts). Therefore r' is in some R(n',l — 1). The weight n' can be easily 
calculated from to and u. 

Our decomposition is clearly unique by definition. □ 

As an example, we consider the decomposition of the reduced form r = 
(1,2, 0,1,0,1) in Figure [2] Clearly we have r = (2, (1,0,1), 1), (0,1). 

Since r' in the decomposition is also a reduced form, we can apply the 
decomposition recursively. An example of a full recursive decomposition can be 
found in Figure[3j Different colors are used for different levels of decomposition. 
Graphically, we can see that at each level, we always have a skewed strip of 
thickness to on the left of position l, and some “dust grains” corresponding to 
u on the right of position l. 

2.3 Construction of generating sequences 

Before proceeding to counting and generating SPM(n), we provide a construc¬ 
tion of a sequence of applications of the FALL rule that allows to obtain any 
given accessible SPM configuration from the initial one. This construction is 
essentially an alternative proof of one direction of Theorem 12.II 
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r = (4,4, 3,2, 2,3, 3,2, 3, 0,0,1,1) 

Figure 3: An example of full recursive decomposition of reduced form 

Definition 2.4. (Generating Sequence) For t £ SPM(n), we say that a finite 
sequence (a*) = ao, ai,..., aj_i generates t if t can be obtained from (n, 0,...) 
by successively applying FALL at index a, for i from 0 to l — 1, which implies 
in particular that each such application must be valid. In this case, we call (a*) 
a generating sequence of t. 

We can say that a generating sequence (ai) is a certificate that t is in SPM(n), 
as it provides a path to go from (n, 0,...) to t by applying FALL. Given a finite 
sequence (ai), we can check its validity by applying FALL accordingly, and if 
it is indeed valid, we also obtain the corresponding t. Conversely, we will now 
provide a method that, given t € SPM(n), constructs a generating sequence 
(a^ of t. We also denote by a the m-folcl repetition of a. We recall that the 
concatenation operation of two sequences is denoted by • 

First we will construct a generating sequence for staircases s(k). Let a, be 
the sequence 0,1,...,* — 1 and be the sequence a* • a»_i •... • a\ (here /cdot 
denotes the concatenation of sequences). Let s(k) be the partition obtained 
by adding k + 1 grains to the first column of s(k). It can be verified that 
Pk applied to s(k) produces s(k + 1). With this observation, we see clearly 
that /3i • /?2 • • • Pk -1 generates s(k) when we start from the initial configuration 
(n, 0,0,...) with n = k(k + l)/2. For n > k(k + l)/2, the same sequence gives 
the configuration ( n — k(k — l)/2, k — 1 , k — 2, ..., 1 , 0), corresponding to the 
reduced form (n — k(k + l)/2, 0,..., 0). 

We want to study the effects of the rule FALL on reduced forms. We define 
the rule FALL' operating on reduced forms r £ R(n, w) as follows: 

r = (r 0 ,.. .,r w ) -» r' = (r 0 , ■ ■ .,n - l,n+i + 1 ,.. .,r w ), 
if n > r i+ 1 + 1. 

This is equivalent to say that the following diagram commutes: 
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FALL (l) 


red w 


To construct t £ SPM(n), we denote by w = sw(s) its staircase width and 
r = red w (t) its reduced form. First we use /?i • 02 • • • • • /3 w -i to construct 
the socle of t, then we pick up the viewpoint of reduced form to construct the 
rest. The remaining task is then to construct a path from the reduced form 
(n — w(w + l)/2,0,..., 0) to r using FALL'. 

For simplicity, we denote by (0) the only element in R{ 0,0). Using Theorem 
El we now define recursively a function Path n<w such that, given a reduced 
form r £ R(n, w), constructs a path from (n — w(w + l)/2, 0,..., 0) to r in the 
Hasse diagram of SPM(n) by applying FALL 7 . 

For w = 0, Path ntW (r) = (), the empty sequence. For w > 0, let r = 
((t, to), l, u ) be the decomposition of r. For i\ > ... > if. > l such that = 

1 for 1 < j < k , we define seqo{u) = a^ ■ at 2 ■ ■ a.i k . Each a tj sends a grain 

( i.e. adds 1) to the component r tj = = 1. In order to construct the to 

layers of grains from position 0 to position l — 1, we only need to repeat m times 
the sequence cq_i • cq _2 ... a.\. We define seq±(l, to) = (ai -1 • cq _2 • ... • a\) m . 
Finally we define Path n:W (r) recursively as follows: 

Path n , }W (r) = seqo(u) ■ seqi(l,m) ■ Path n -k-im,l-i(t)■ 

We have the following proposition stating the correctness of the construction 
P ath HiW . 

Proposition 2.7. For s £ SPM(n ) with w = sw(s), [i\ ■ $2 ■ ■ /3 w -i ■ 

Path niW (red u ,(s)) is a generating sequence of s. 

Proof. We perform an induction on w. The base case w = 0 is trivial. To pro¬ 
ceed by induction, we suppose that the proposition is true for any s £ SPM(n) 
having width smaller than w. By combining the notation in this proposition 
and the result of Theorem 12.21 we have red w {s) = (l, u, m),r'. We observe that 
r' is a reduced form of a certain configuration s' with width l < w. It is easy to 
verify that the sequence Path UtW (by its recursive definition and our induction 
hypothesis) constructs the correct reduced form, which completes construction 
of the socle obtained by the successive applications of the /Vs. □ 

This proposition can be seen as a constructive proof of one direction of 
Theorem I2TT1 We should be aware that this construction is not unique. For 
instance, in the recursive definition of Path niW , by exchanging the order of 
seqo(u) and seqi(l, to), we can get a different valid construction. From a physical 
point of view, it may be interesting to study the number of generating sequences 
of configurations in order to have better understanding over the generic evolution 
of an SPM model. 







2.4 Recursive formula for |SPM(n)| 

The recursive structure of SPM(n) described by Theorem 12.21 can be used to 
give a counting formula for SPM(n). 

We define c(p, w) = \R(p+w(w+l)/2,w)\. The following proposition follows 
directly from the definitions of c{jp,w) and R(p,w). 

Proposition 2.8. For a natural number n, we have 



W>1 

w(w-\-l)<2n 


The reason we choose c(p,w) = \R(p + w(w + l)/2,w;)| is that for all r £ 
R(j> + w(w + l)/2, w), we have ]TN r* = p. Here p represents the number of 
grains located “above” the socle s(w) = (w, w — 1,... ,1,0). As a consequence 
of Theorem 12.21 we have the following recurrence for c(p, w). 

Theorem 2.3. The value of c{p,w) is uniquely determined by the following 
recurrence. For w > 0, we have c(0,w) = 1. For p ^ 0, we have c(p, 0) = 0. 
For the remaining cases, 



Proof. This recurrence comes directly from the decomposition of reduced forms 
r = (( t,m),l,u ). The base cases of the recurrence can be easily verified. The 
summation index l (resp. m) corresponds to the integer l (resp. m) in the 
decomposition. The summation index i stands for the number of parts of u 
equal to 1. Binomial coefficients arise by taking the number of all possible 
sequences u of 0’s and l’s having exactly i parts equal to 1 (and all other parts 
equal to 0, as in Theorem E2). The special case l = 0 is treated in the first 
term on the right hand side. This is the case where t is empty. The p grains in 
r can be only placed on w columns in u, and at most one grain can be placed 
on each column. Therefore, the first term is non zero only when 0 < p < w. □ 

Now we want to evaluate the complexity of computing |SPM(n)| using this 
recursive formula. 

All binomial coefficients (f) needed to calculate c(n,w) verify a < w. By 
Proposition CO we have w < \[2n. Therefore, by memorizing results (in the 
manner of dynamic programming), we can precalculate all of them using Oin ) 
additions. 

In the recurrence for c{p,w), we notice that m cannot exceed p/l , therefore 
for a fixed l < w, there are at most p/l possibilities for the value of m, thus 
we have Y^i=iP/^ = ?Ei=i = 0(plog(w)) possible pairs ( l,m ). Since 1 < 
i <w, in the recurrence for c(p,w) there are at most 0(wp log(u)) terms, thus 
0(wplog(w)) arithmetic operations are needed to calculate each c(p,w), given 
the value of all c{p ', w') with p' < p and w' < w. The total number of arithmetic 
operations for calculating all c(p,w), for all p < n, is bounded by: 




p<n 1 <w<V2n 
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It follows from Proposition 12.81 that we only need 0(n 3 log(n)) arithmetic 
operations to compute |SPM(n)|. According to Tj, we can bound |SPM(ra)| by 
c n (for a certain constant c), thus all coefficients involved have 0(n) bits, thus 
we know that we can compute |SPM(?r)| in 0(n 4 log 2 n log log n) time using fast 
integer multiplication. 

As a remark, given the recursive formula for c(p, w ), it is straightforward to 
construct a uniform random generator of SPM(n) by computing all c(p, w) and 
generate configurations recursively in a uniformly random way using appropriate 
probabilities computed with c(p, w). 

2.5 A CAT algorithm for SPM(n) 

It is clear that the exhaustive generation of SPM(n) reduces to the exhaustive 
generation of SPM(n, w) with the staircase width w varying from 1 to [\/2nJ, 
which in turns reduces to the exhaustive generation of reduced forms in R(n , w). 
The unique decomposition of reduced forms in Theorem ^. 21 thus gives a natural 
way to exhaustively generate reduced forms in a recursive fashion. Essentially 
our algorithm will be an algorithmic transcription of Proposition 12. (SI and of 
Theorem 12.31 Algorithm [T] is an example of such a transcription. It should 
be called initially with d = 0. The cases on l are for the further complexity 
analysis. 

We now explain the data structure we use. From the notation in Theo¬ 
rem E1 the unique decomposition of reduced form gives each reduced form r 
an expression as a list of triplets ( lo , uq, m o),..., (Id, Ud, nid) with natural num¬ 
bers k > 0, nii > 0 and (0, l)-sequences zq, with the condition that the sequence 
(^)o<i<d is strictly decreasing. We will adopt this notation for our exhaustive 
generation algorithm. In the generation, the natural numbers li,rrii come natu¬ 
rally from loop indices, and the rest consists of the generation of zq, which are 
(0, l)-sequences, with given length and given weight (be. total number of Is). 
There are various CAT algorithms for generating (0, l)-sequences with given 
length and weight, for example those presented by Knuth in Chapter 7.2.1.3 of 
[U and by Ruskey in [8] in Section 4.3. Any one of these methods can be used 
as a subroutine to generate u. Many such algorithms have a linear initialization 
when the given length is equal to the given weight, and they might fail to be 
CAT in this special case. However, this can be fixed by using a boolean variable 
associated to zq to indicate this special case. Therefore, we can consider all 
operations in Algorithm |T] to be performed in constant time. 

We now analyse the time complexity of Algorithm |T] According to Section 
4.3 in (HI , we only need to analyse the form of the recursion tree. For each call 
of Generate(p,w,d), it is clear that at least one of the generated SPM config¬ 
urations is produced immediately in this call, by putting all the p grains into 
the first column. Therefore, the number of nodes of Generate(p,w,d) in the re¬ 
cursion tree of Generate(n — w(w — l)/2,u; 1 0) is bounded by |SPM(n, w)\, thus 
the total number of nodes is bounded by 2|SPM(n, w)\. We also know that, in 
Algorithm [TJ the time spent to spawn a child for each node in the recursion tree 
is bounded by a constant. It is then immediate that our algorithm is CAT for 
SPM(n,w), thus also CAT for SPM(n). 

We now analyse the space complexity of Algorithm [1] first expressed in 
terms of memory cells, then in terms of bits. We notice that generating u of 
length k uses 0(k) extra memory. As the total length of all u in a recursive 
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i Generate^, w , d) 

Result: Each call yield() (including those in recursive calls) sees a new 
r £ R(p + w(w + l)/2, w) in the array A of triplets ( li , u,;, m*). 
Throughout the code, l is the position of first zero, i the number 
of Is in u and m the minimum before the first zero. 
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begin 

if p = 0 then 
L yieldQ; return; 

// Case l = 0 

if p < w then 

foreach u of length w with p Is do 
|_ A[d\ <- (0,0, u); yield(); 

// Case l = 1 

for i 4 — 0 to min(u; — 1 do 

foreach u of length w — 1 with i Is do 

L A i d ] ^ *); y ield 0; 
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// Case l > 2 
for Z 2 to w do 

for i i — 0 to min(w; — /,p — I s ) do 

foreach u of length w — l with i Is do 

for mi— 1 to L^J do 

^ A[d] i— (l, it, m); Generate^ — * — Im, l — 1, d + 1); 
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return; 


Algorithm 1: Recursive Generation of R(p + w(w + l)/2 ,w) 


call is bounded by w, we know that this part of memory consumption is 0(w). 
Secondly, w decreases at each recursive call, therefore the recursion depth is 
0(w). Since Generate(p,w,d) only uses a constant number of scalar variables 
besides the array A, the total stack memory consumption is 0(w). Adding 
the memory needed for the array A and for generating it, the total number of 
memory cells used in Generate(p,w,d) is 0{w). As the value of each memory 
cell is bounded by max(n,w), space complexity of generating SPM(n, w) is 
0(w log(max(n, w))) bits. For the generation of SPM(n), a simple reuse gives a 
total space complexity of Oi^/Ti log(ri)) bits. 

3 Generalization to ice pile model 

We will now generalize previous results to the Ice Pile Models IPMfc(n), using 
the same terminology as in previous sections. 

For IPMfc(n), we define the following staircases for w > 0 and 1 < l < k: 

s(w,l) = ( w,... ,w , w - 1 1) . 

I k k 

In Figure^] several examples of staircases for k = 2 are presented. 
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*(3,2) 


*(4,1) 


□ m IThn ITT! 

a(l,l)a(l,2) s(2,1) s(2,2) s( 3,1) 

Figure 4: The first staircase bases of IPM 2 


These staircases are clearly stable by all rules of IPM* (n). We define anal¬ 
ogously the staircase basis Bk = {s(w,l) \ w > 0,1 < l < k}. It is clear that 
s(w, l ) < s(w', l') (sequence order) if and only if (w, l ) is not larger than (w' , V) 
in lexicographical order. We can see that < is a linear order over Bk- For 
t G IPMfc(n), we define sw(t) = ( w,l ) such that s(w,l ) < t, and for any s G B 
with s < t, we have s < s{w,l). For example, for t = (8, 8,5,5) G IPM 2 (26), 
we have sw(t) = (2, 2) and s(sw(t)) = (2,2,1,1), while, if the same t is seen as 
a configuration of IPM 5 (26), we have sw(t) = (1,4) and s(sw(t)) = (1,1,1,1). 
As before, with respect to the sequence order, s(sw(t)) is the largest staircase 
among all those that are smaller than t. It is clear that we cannot apply SLIDE*, 
on any s(w,l). We also have the following analogue of Proposition 12.21 

Theorem 3.1. For t,t' G IPMk{n ) and t —> t', we have s(sw(t)) < s(sw(t')). 

Proof. Set b = s(sw(t)). By definition of sw, it suffices to prove that b < t'. 
Suppose that t' is obtained from t by applying the rule SLIDE* on column c. 
Therefore, t c = t' c + 1; for some suitable p < k, we have t c + p = t' c+p — 1 and 
t c +p <t c — 2; and for all j $ {c, c + p}, we have t 3 = f'. By definition, b < t. 
The only column that may prevent b < t' is column c. However, b c+p > b c — 1 
for any p < k (by definition of 5), thus t c + P > b c+p > b c — 1. Since we also have 
t' c = t c — 1 > t c+p + 1, we have t' c >b c , which concludes the proof. □ 

We now propose a few definitions similar to those settled for SPM. 

Definition 3.1. Staircase width and reduced form for IPM 

We define IPM*(n, w, l) as the subset of IPM*(n) of all elements with stair¬ 
case width (w,l), or formally IPM*(n,u;, l) = {s G IPM*(n)|sw(s) = (uqZ)}. 
We can see that the family {IPM*(n, w, l)} w ,i is a partition of the set IPM*(n). 

For s G IPM*(?z, w, l), we say that red( w j)(s) = (s, — b(sw(s))i)i> 0 is its re¬ 
duced form. By definition, every reduced form is a sequence of natural numbers. 
We denote by f?*(n, to, l) the set of reduced forms of elements in IPM*(n, to, l). 

Some examples of IPM configurations and their reduced forms are illustrated 
in Figure [21 

To obtain analogues of the decomposition theorem, we start from the char¬ 
acterization of elements in IPM*(n). The following characterization of IPM(n) 
is first given in [3], then in .7] it is used to give an exhaustive generation al¬ 
gorithm for IPM(?t). We adapt the following notations from [7|. We denote 
by 1 the sequence (p,...,p) of n elements equals to p and we recall that the 
concatenation operation of two sequences is denoted by ■. 

Theorem 3.2. A partition s is in IPMk(n) for a certain n if and only if it does 
not contain the following forbidden patterns (for p > 0 and h > 1): 

• p[k ~*" 2 1 
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s 


5 


= (7, 7,4,4,0,...) 
Basis: s(2,2) 
red {2 , 2 )(s) = (5,5, 3,3,0,...) 


red (4j i)(s) 


( 6 , 4 , 4 , 3 , 2 , 1 , 1 , 1 , 0 ,...) 

Basis: s(4,1) 

= ( 2 , 1 , 1 , 1 , 0 , 0 , 0 , 1 , 0 ,...) 


Figure 5: Some examples of reduced form in IPM 2 


• (p+ l)[ fe+1 l • 1 

• (jp + ft.)I fc+1 l • Uti(P + h — i)I fc l • p[ fc+1 l 

Using this characterization, we will prove an analogue of Proposition ^. 51 for 
the ice pile model. 

Lemma 3.1. The largest (in the sequence order <) IPM k configuration s with 
Sq < w is w ■ s(w , k). 

Proof. Let s be an IPM^. configuration with so < w. We have s, < w for all i > 
0. Theorem 12. II implies that, for any integers i > 0 and p > 0, s^+pfc+i < Si — p, 
thus for any 0 < i < k and p > 0, s p k+i+ 1 < w — p. We observe that this is 
equivalent to s < w ■ s(w). We conclude by noticing that w ■ s(w, k ) is also an 
IPMfc configuration. □ 

Proposition 3.1. For s € IPMk(n,w,l), we have s; + fc( w _i) + i = 0. 

Proof. For s £ IPMfc(n, w, l), we have sw(s ) = (w,l) and we denote by b = 
s(sw(s)) the staircase base of configuration s. By definition of b , there is i = 
l+kp > 0 for a certain p such that s, = bi = w—p—1. Since suffix (s,, Sj+i,...) is 
also an IPM^. configuration, by applying Lemma 13.II we know that it is smaller 
than w p 1 ■ s(w - p - 1, k), thus sj+fc( tt _i)+i = s i+(u ,_ p _i) fc+ i =0. □ 

This proposition means that every reduced form r £ Rk{n,w,l) is in fact a 
(l + k{w — 1) + l)-tuple of natural numbers. We will now characterize elements 
in Rk{n,w,l) by the following analogue of Proposition 12.61 

Proposition 3.2. A (l + k(w—l) + l)-tupler = (ro,ri,..., U+fc(w-:L)) of natural 
numbers is in Rk{n, w, l) if and only if 

• fi) r i = n - Iw - kw(w - l)/2; 

• (ii) There exists an index 0 < io < l + k(w — 1) of the form l + kpo (for 

some integer po) such that ri 0 = 0; 

• (Hi) For all i > 0 of the form l + kp (for a certain integer p), we have 
U > r i+ 1 > ... > r i+k -1 > r i+k - 1; 

• (iv) For all i > 0 and j = i + kp + 1 (for a certain integer p > 0), we 

have Vi > rj — 1 when i = l — 1 (mod k), and r.j > rj otherwise. 
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Proof. Let s £ IPMfc(n, w, l) and r = redt w ^{s). Conditions (i), (ii) and (iii) 
come from the definitions of IPMfc(n, w , l) and of redr Wt iy For the condition (iv), 
we notice that Theorem ld. 21 implies that for any integers i > 0 arid j = i + kp+l 
(for a certain integer p > 0), Sj < Si —p. Since we have s(w, l)i = s(w , l)j +p +1 
when i = l — 1 (mod k) and s(w,l)i = s(w,l)j + p otherwise, we easily verify 
that condition (iv) holds. 

Conversely, let rbea(J + k(w — 1) + l)-tuple that verifies conditions (i), (ii), 
(iii) and (iv), and t = r + s(w, l). It follows from conditions (i) and (iii) that t 
is a partition of n, and it suffices to prove that t £ IPM& (n) , because condition 
(ii) will ensure that t has the correct basis. Suppose that there is a forbidden 
pattern in t between column c and c + kp + 1 for some p > 0, and we have 
t c+k p+i = t c -p+l. If c = l — l (mod k), we have s(w,l) c = s(w,l) c+kp +i+p+l, 
thus r c = r c +fcp+i — 2; otherwise, we have similarly r c = r c +kp +1 —1. This cannot 
happen when (iv) is verified, thus t £ IPMfc(n). □ 

We notice that conditions (i) and (ii) ensure that the reduced form is in 
the Rk(n,w,l) with correct parameters. Conditions (iii) and (iv) are stable by 
prefix-taking, the same as in the case of SPM. However, these two conditions, 
and also the length of tuple, are parametrized by w and l. Simply taking prefix 
will preserve (iii) and (iv), but with parameters w, l not compatible with the 
length of tuple. Therefore, if we mimic the decomposition theorem for SPM in 
a naive way, the part before the first zero will not be a valid reduced form. We 
try to circumvent this problem by extending our definition of reduced form. 

Definition 3.2. Extended reduced forms, augmented reduced forms 

For a pair of positive integers (w, l ) and a tuple t of length l + k(w — 1) + 1, we 
say that t is an extended reduced form if t verifies conditions (ii), (iii) and (iv) 
in Proposition [321 and t is called an augmented reduced form if conditions (iii) 
and (iv) are verified. 

We denote by R' k (w,l) the set of extended reduced forms, and Ak{w,l ) the 
set of augmented reduced forms. Clearly we have R' k (w,l) C A k {w,l). 

Clearly the subset of all (l + k{w — 1) + l)-tuples in R' k (w , l) is exactly the 
union of Rk{n , w, l) for all possible n. 

We now introduce a function that will be used to turn augmented reduced 
forms into the more regular extended reduced forms. 

Definition 3.3. For positive integers k and l such that 0 < l < k, we define the 
function pl l k on tuples with arbitrary length of non-negative integers as follows: 
set r' = pl l k {r) : we define r' = r, — 1 for i = l (mod k) and r' = otherwise. 
This function is undefined when there exists some i = l (mod k) such that 
n = 0 . 

The example in Figure [S] shows graphically the effect of pl l k . Intuitively, if 
r is the augmented reduced form of some IPM^ configuration with respect to 
the basis s(w,l), then pl l k {r ) is the augmented reduced form of the same IPM^ 
configuration with respect to the next basis, in the linear order for B k . 

We now investigate some properties of pl l k in the following lemma. 

Lemma 3.2. The function pl l k verifies the following properties. 

1. The function pl l k is undefined on the set of extended reduced form R' k (w , l). 
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r = (4,3, 3,2, 0,1,1,0,1,1) € A 3 (3, 2) r' = (4, 3,2, 2, 0, 0,1, 0, 0,1) G A 3 ( 3,3) 
Figure 6: An example of application of pl l k , with k = 3, l = 2 


2. For a tuple r G A k (w,l), if pl l k (r) is defined, we have pl l k (r ) G A k {w' ,1'), 
where w' = w + 1, l' = 1 if l = k, and w' = w, V = 1 + 1 otherwise. 

3. For a tuple r G A k (w,l), we recursively define the sequence of tuples r = 

r (°), 7 ~(i),... by r( l+1 ) = pl k ,] (rW) G A k (w^ +1 \l^ +1 ^). This sequence 
becomes undefined after a certain index c satisfying r^ G R k (w^ c \l^). 
Moreover, does not depend on r. 

4- In the case of the previous assertion, we say that r = Aug Wi i(r^ c \c) is 
equal to r ^ augmented by c. Regarded as a function, Aug Wi i is a bijection 
between {( rc)\r' G R k (w^ ,1^), c G N} and A™’ 1 . 

Proof. The first assertion follows from the definitions of both R k (w,l ) and pl l k . 
More precisely, the elements of R’ k (w,l) verify condition (ii) in Proposition [321 
which prevents pl l k to be defined. 

The second assertion comes from simple verification of conditions (iii) and 
(iv) in Proposition 13.21 for pl l k (r). 

We now prove the third assertion using the result of the second. We will 
first prove that the sequence terminates, then discuss the properties of r^ and 
( w^i ), ?W). 

To show that the process terminates, we notice that pl l k is always a decreasing 
function with respect to the sequence order, and strictly decreasing in the case 
l = 1, since a tuple must have its first element. If the sequence r^ = r, ... 
does not terminate, the case l = 1 will occur an infinite number of times, thus we 
can extract an infinite strictly decreasing sequence from the original one. This 
contradicts the well-foundedness of the sequence order of tuples of non-negative 
integers with fixed length. Thus termination of the process follows. 

For the iterative process to terminate, pl l k must be undefined on A c \ that 
is to say condition (ii) is verified, following the same reasoning as in the first 
assertion. Combining with r^ G A k (w^ c \l^), we have r ^ G R' k (w^°\l^). 
The independence of (w^\l^) from r is implied by the independence of w',V 
from r in the second assertion. 

For the last assertion, Aug Wt i is clearly surjective. We also notice that, given 
a pair (r', c) with r' G R'AvM’,1^), it is easy to uniquely reconstruct a tuple r 
such that r = Aug Wt i{r^ c \c) by reversing the recursive process indexed by the 
sequence (w^°\ 1^,1^) independent of r'. Therefore, A k ' 1 is also 
injective, which proves the assertion. □ 

We notice that s(w,l ) < s(w',l') (as defined in the second assertion of 
Lemma l3.2ll are consecutive elements in the linear order B k with respect to the 
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sequence order. In fact, let m be the length of tuple r, and s(w, l)\ m the prefix of 
s(w, l) of length m, then we can easily verify that r + s(w, l)\ m = r' + s(w', V )| m , 
where addition is intended as pointwise. This equality means that the function 
pl l k transforms an augmented reduced form on a certain socle to another aug¬ 
mented reduced form on the smallest socle that covers the previous one, and 
these two augmented reduced forms are equivalent in the sense that they actu¬ 
ally give the same prefix of a configuration, but with the removal of different 
staircases. 

With all these modifications, we can state an analogue of our SPM reduced 
form decomposition theorem. 

Theorem 3.3. An extended reduced form r £ R k (w,l) can be uniquely decom¬ 
posed into the following form: 

T — (^0j * * • j tl-\-)zp— 1,0, Uq, • * - 5 ^m) j 

with some integer p such that t = (to ,..., ti + k p ~ i) and u = (uq, ..., u m ) verify 
the following conditions: 

1. If u is not empty, we have Ui £ {0,1} for i = — 1 (mod k), and Ui = 0 
otherwise. 

2. We have t £ Ak(w, l), but t ^ R'j-(w, £)■ 

This decomposition will be denoted as r = ((t 1 ,c),p,u), with t = Aug w g(t^ c \c) 
and t! = t^ c \ When t is empty, we take c = 0. 

Proof. The existence of index i = l + kp such that r* = 0 is given by condition 
(ii). To ensure t R’ k (w , l), we take the smallest such index. By condition (iii), 
we know that iq £ {0,1} for i = —1 (mod k) and iq = 0 otherwise, and this 
does not violate condition (iv). It is clear that t £ Ak(w,l), since condition (iii) 
and (iv) are invariant under prefix-taking. Therefore, from the third assertion 
of Lemma [3.21 we have the existence of (t^ c \c) as a representation of t. We 
now have the validity of this decomposition. The uniqueness is provided by the 
uniqueness of smallest index i = l + kp such that r* = 0, and the uniqueness 
of the pair (r^ c \c) giving r = Aug Wi i(A c \c) provided in the last assertion of 
Lemma 13.21 □ 

Since Rk(n , w, l) is the subset of R' k (w, l) consisting of all (,l + k(w — 1) + 1)- 
tuples with correct weight, we can use this decomposition theorem of extended 
reduced forms to enumerate and to generate IPM configurations, following the 
same approach as previously done for SPM. 


4 Future work 

We would like to extend this approach to more general sand pile models, for 
example BSPM, a bi-dimensional version of SPM. However, this extension does 
not seem to be easy. There are several difficulties. Firstly, rules now involve 
two directions, which weakens the foundation of staircase bases on well-behaving 
rules. Secondly, we do not yet have a good characterization of configurations in 
BSPM, even for stable ones. Lastly, simulations show that fixed points in BSPM 
have a great variety of different shapes, which would be difficult to approximate 
using a “small” set of staircase bases. 
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